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ABSTRACT 


^Unsteady  flow  about  the  well-rounded  nose 
of  a subsonic  airfoil  is  investigated  from  the 
viewpoint  of  leading  edge  separation.  For  an 
airfoil  undergoing  forced  pitching  and  heaving 
motions  in  a uniform  flow,  HhaP  fluid 
accelerations  about  the  leading  edge  can  be 
enormous  According  to  inviscid  flow  theory. 
Such  accelerations  are  limited  by  viscous  flow 
and  separation  realities. 


The  method  of  matched  asymptotic  expansions 
r-.  is  used  to  develop  a uniformly  valid  first  order 
approximation  to  the  inviscid  flow  about  the  — 
— ijctetV  -The-  invi  scid " flow  aootrt— the^airfoTl'T 
jVvtk/Teading  edge. is  driven  by  a hi  story -dependent 

term  related  to  the  airfoil's  transverse 

motions.  Applying  this  flow  to  the  laminar 
boundary  layer  flow  at  the  airfoil  nose  produces 
possibilities  for  a laminar  boundary  layer  to 
separate.  * nuT  ,./m  , 

f The  Nbore-Kott- Sears  ( H-R-S)  conditions  for 
unsteady  boundary  layer  separation  do  not  appear 
to  be  useful  for  this  problem.  A methodology  is 
''^proposed  for  predicting  leading  edge  dynamic 
stall  based  upon  relating  properties  of  the 
envelope  of  the  unsteady  part  of  the  boundary 
layer  speed  and  shear  stress  to  the  steady 
part  of  the  boundary  layer  flow.  ^ The 
development  is  proposed  as  a tool  \ for 
determining  the  useful  limit  for  applying 
attached  inviscid  airfoil  flow  theory. 
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NOMENCLATURE 

complex  amplitude 
acceleration  vector 
suction  strength 
Fourier  velocity  coefficient 
constant  related  to  on-set  flow 

■ Re  0,  Theodorsen 
« lm  0,  Theodorsen 
complex  velocity  potential 
Bernoulli  constant 
displacement  function 
unit  complex  number,  space 
unit  complex  number,  time 
steady  viscous  function 
modi  fled  Bessel  function 
local  pressure 

flow  spend 

Reynolds  number  based  on  chord 
nose  radius  of  curvature 
time 

on-set  speed 
tangential  speed 
velocity  components 
velocity  components 
coordinate  pair 
coordinate  pair 

■ X + 1Y,  complex  position 
coordinate  pair 

■ V 1c.,  complex  flow  function 

■ v^T  y,  similarity  parameter 


= ait 

Theodorsen  function 
= 8 + ia,  complex  variable 
heaving  parameter 
pitching  parameters 
fluid  density 

reduced  frequency  of  the  motion 
acceleration  potential 
velocity  potential 
stream  function 
kinematic  viscosity 
circular  frequency 


INTRODUCTION 

Based  upon  the  linear  aerodynamic  theory, 
Wu*  developed  hydrodynamic  analyses  for  optimum 
pitching- heaving  motion  of  a rigid  wing.  The 

optimum  problem  was  to  minimize  the  time 
averaged  energy  loss  coefficient  Cg  under  the 
constraint  that  the  time  averaged  thrust 

coefficient  CT  was  fixed.  That  is,  for 

C_  ■ "C-  >0,  what  phasing  between  the 

puchinjSheaviny  motions  minimizes  the  rate  at 
which  energy  is  lost  to  the  flow  in  the  shedding 
of  kinetic  energy  at  the  airfoil's  trailing 
edge?  The  time  averayed  thrust  coefficient  is 
comprised  of  a mean  thrust  delivered  from  the 
plate  surface  and  a mean  suction  due  to  the 
inviscid  flow  accelerating  about  the  airfoil 
leadlny  edye.  Uu  determined  that  the  ratio  of 
mean  suction  thrust  coefficient  to  total  thrust 
coefficient  Cc/Cj  has  a minimum  at  a reduced 
frequency  of  the  motion  o * o^CC-  ) • Outside 
of  the  region  of  o » o , the  suction  force  can 
become  so  large  thatm  leadlny  edge  stall  Is 
Inevitable.  The  present  paper  advances  a 
methodology  for  determining  the  limit  on  the 
attainable  Cc  for  oscillatory  motion  of  a thin 
airfoil  based  upon  the  behavior  of  a laminar 
boundary  layer  near  the  airfoil's  leadlny  edge 
stagnation  point. 

Sychev's  impressive  manuscript^  shows  that 
under  certain  restrictions  on  the  acceleration 
of  the  flow,  the  point  of  separation  Is  Inviscid 
in  nature.  This  result  contrasts  with  the 
steady  separation  problem.  Furthermore,  the 
unsteady  separation  point  is  not  situated  on  the 
surface  of  the  body. 

The  approach  taken  herein  has  been  to  treat 
the  leadlny  edye  separation  problem  in  the 
context  of  Its  relationship  to  the  full  flow 
about  the  oscillating  airfoil.  In  treating  the 
problem  In  this  way,  attention  is  given  to  the 
situation  that  leading  edge  separation  appears 
to  be  a local  phenomenon  in  the  sense  that  it 
occurs  in  a region  where  the  airfoil  has  large 
curvature,  and  experiences  a stagnation  point 
flow.  Once  dynamic  separation  has  occurred,  the 
aerodynamic  theory  used  to  describe  the  attached 
flow  about  the  airfoil  no  longer  applies. 


The  approach  taken  connects  the  linear 
airfoil  theory  and  the  unsteady  oscillating 
laminar  boundary  layer  theory.  The  latter  has 

been  developed,  for  example,  in  References  [3]  - 
[7]. 


It  is  well  known  that  the  linear  airfoil 
theory  breaks  down  at  the  airfoil  leading 
edge.  This  deficiency  has  been  eliminated  in 
this  development  by  asymptotic  matching  with  a 
local  inviscid  unsteady  solution  about  an 
osculating  parabola.  Results  from  matching  then 
open  the  way  for  investigating  the  behavior  of 
the  laminar  boundary  layer  in  the  region  of  the 
airfoil  nose.  This  boundary  layer  is  driven  by 
the  gross  oscillatory  motion  of  the  airfoil  and 
by  its  steady  forward  speed.  The  unsteady 
driver  turns  out  to  be  directly  related  to  the 
suction  strength  about  the  airfoil's  leading 
edge,  (as  derived  from  the  inviscid  theory). 


v(x,t) 


-5T+  U15T 


]h( x ,t) 


< 1 


v t > o 

For  heaving  and  pitching  motion  of  a rigid  plate 
at  the  frequency  w,  the  displacement  function 
h(x,t)  can  be  expressed  as 


h(x,t)  = [50/2  + (51  + j52)]exp  (jut) 

where  |x|  < 1 and  £ , 5,,  and  ^ are  real.  The 
pitching  axis  is  at  the  midchord,  x = 0.  The 
amplitudes  of  pitching  and  heaving  are 
k +js  | and  5 /2,  respectively.  Pitching  leads 
the  . heaving  0 motion  by  a phase  angle 
tan"1^/^). 

A convenient  form  of  expressing  v(x,t)  is 
by  its  Fourier  cosine  series.  This  series 
contains  only  two  two  terms  for  the  specified 
transverse  motion.  That  is. 


Flat  plate  unsteady  boundary  layer  theory 
is  applicable  in  this  Investigation  whenever  the 
boundary  layer  thickness  is  small  compared  with 
the  radius  of  curvature  of  the  airfoil's  nose, 
lie  tacitly  assume  this  to  be  the  case. 

KINEMATICS 

lie  consider  the  small  amplitude  heaving 
and  pitching  motion  of  a thin  symmetric  airfoil 
in  a steady  uniform  stream  U_.  To  describe  the 
kinematics  of  such  motion  we  introduce  a 
Cartesian  coordinate  system  (x,y)  with  the  x- 
axis  aligned  with  the  airfoil's  mean  chord 
line.  The  direction  of  the  free  stream  is  along 
the  positive  x-axis.  The  airfoil's  transverse 
displacements  then  occur  along  the  y-axis  and 
are  prescribed  by  a function  h(x,t)  of  chordwise 
position  x and  time  t.  Figure  1 Illustrates  the 
transverse  displacement  of  a typical  wing 
section  with  respect  to  the  (x,y)  coordinate 
system. 

The  inviscid  wing  boundary  condition 
requires  the  normal  velocity  of  the  wing 
relative  to  the  (x,y)  coordinate  system  be  equal 
to  the  normal  velocity  ^»n  of  the  fluid  adjacent 
to  the  wing.  Here  _n  is  the  unit  outward  normal 
vector  on  the  wing  and  ^ is  the  fluid  velocity 
relative  to  the  inertial  reference  frame 
resolved  into  components  (U+u,v)  along  the  (x,y) 
axes. 

Neglecting  products  of  small  quantities 
compared  with  those  occurring  linearly,  the 
kinematic  boundary  condition  specifies  the  v- 
coraponent  of  the  fluid  velocity  adjacent  to  the 
wing.  Consistent  with  approximations  already 
made,  this  component  can  be  specified  along  the 
x-axis,  giving 


Airfoil 


v(x,t)  * bQ(t)/2  + b^(t)  cos  e 

(1) 

b0(t)  = uf2Cx+  ]eJut 

(2) 

bjU)  s -Uo(s2  - j^le^ 

(3) 

where  x » cos  0,  9 * iot  and  o =■  u/u 

is  the 

reduced  frequency  of  the  motion  based  on  the 
unit  half-chord. 

INVISCID  DYNAMICS 

In  an  incompressible  flow  field  devoid  of 
external  forces  and  Internal  viscosity,  the 
principle  of  conservation  of  mass  leads  to  the 
expression  * U.  Conservation  of  rectilinear 
momentum  leads  to  the  Euler  equation  herein  the 
pressure  gradient  is  balanced  by  the  fluid 
acceleration  a * -p'l£p.  This  equation  is  valid 
In  any  Inert fal  reference  frame.  The  absolute 
acceleration  a.  measures  the  rate  of  change  of 
following  a particle,  p is  the  fluid  density, 
and  p Is  the  local  Instantaneous  pressure.  If 
in  addition,  the  flow  field  is  irrotational 
then  v x ^ * U and  a scalar  function  * exists 
such  that  ^ * U+V*  and  « 0.  Here  the  flow 
field  is  defin¥d— to  be  the  region  exterior  to 
the  wing  and  its  shed  vortex  sheet. 

An  integral  of  the  momentum  equation  can  be 
obtained  by  substitutiny  jj  * U+v*  into  the 
acceleration  a and  neglecting  products  of  small 
quantities.  Die  resulting  Integral  becomes 

(-jf  ♦ U^l*(x,y,t)  • »(x,y,t) 

where  *,  the  Prandtl  acceleration  potential , 
measures  the  variation  of  the  pressure  from  the 
static  level,  »(x,y,t)  » [p,-  u(x,y,t)]/p. 
Applying  the  Laplace  operator  r to  the  above 
Integral  results  in  r*  « u as  the  relevant 
field  equation  for  the  unsteady  Inviscid 
incompressible  thin  airfoil  problem. 

INVISCID  OUTER  PROBLEM 

The  appropriate  airfoil  boundary  value 
problem  has  been  solved  in  terms  of  the 
acceleration  potential  by  Uu*1  In  an  elegant 
treatise  on  the  hydrodynamics  of  swinwing 
propulsion.  In  particular,  he  has  determined 
that  the  pressure  difference  across  the  wing, 

| x | < 1 is 
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4p  * p*(x,t)  - p+(x,t)  * 2p*+( x ,t ) (4) 

where  4+(x,t)  is  the  acceleration  potential 
evaluated  along  the  x-axis,  y « 0+. 


For  time  harmonic  motion  that  has  persisted 
indefinitely,  the  acceleration  potential 
evaluated  on  the  topside  of  the  wing  is 


where 

V 


(5) 

(6) 

(7) 


aQ(t)  = b1(t'/  - 0(o)[bQ(t)  + b^t) ] (8) 

In  this  last  expression,  0(0)  is  the  Theodorsen 
function,  which  is  expressible  in  terms  of 
modified  K-type  Bessel  functions,  as 


The  simplest  representation  of  the  inviscid 
flow  about  a well-rounded  nose  of  a thin  airfoil 
that  preserves  the  essential  character  of  the 
problem  is  the  flow  about  an  infinite 
parabola.  Upon  magnifying  the  detail  at  the 

leading  edge  what  appears  is  the  flow  about  a 
so-called  osculating  parabola.  Such  a flow  is 
comprised  of  an  on-set  stagnation  point  flow  and 
a tangential  or  parallell  flow.  Figure  2 
Illustrates  the  inviscid  flow  about  a parabolic 
cylinder.  The  geometry  is  described  by  a 

coordinate  system  (X,Y)  with  origin  at  the  oase 
of  the  parabola  and  with  /.  being  the  axis  of  the 
parabola. 


If  we  employ  the  conformal  transformation 
Z = X -x2  (13) 

where 

Z = X + iY  (14) 

x * 6 + ia  (15) 


O(o)  * Kj( ja)/[KQ( jo)  + KjUcr)] 

0(  a)  = F( a)  + jG(a) 

Substituting  Eg.  (1)  into  Eg.  (7)  gives 
AU.t)  * cQ(t)  + Ec^(t)  + C2c2(t) 


(8) 


(1U) 


the  flow  field  in  the  physical  Z-plane  is 
tranformed  onto  the  left  half  x-plane,  as 

presented  in  Figure  3.  That  is,  the  conformal 
transformation  takes  the  parabola  onto  a 
straight  line.  The  upper  branch  of  the  parabola 
goes  to  the  positive  imaginary  x-axis.  The 
lower  branch  goes  to  the  negative  x-axis. 


where 


c0(t)  ■ - t6o  - 61  + Ubo1/Z 

) 

Cj(t)  * - [Ubl  + bQ/z] 

} (ID 

c2(t)  - - bj/2 

) 

Here,  the  dot  denotes  time 

dlfferentatlon 

Substituting  Egs.  (5),  (6),  (1U)  and  (11)  into 
Eg.  (4),  gives  the  pressure  jump  across  the 
wing,  | x | < 1.  That  is, 

£ ■ uaoV5t"+  z(v  xc2,V^"  (i2) 


Eguation  (13)  can  be  used  to  provide  the 
inverse  transformation  yielding  x as  a function 
of  Z.  By  selecting  the  negative  branch  so  that 

< « [1  - Vl  * « j/2  (16) 

then  the  parabola  B * U yields  a - « Y. 

In  terms  of  the  complex  velocity  potential 
function  f,  the  stagnation  point  flow  and  the 
parallel  flow  can  be  readily  represented  in  the 
x- plane  by 

f * - Usx2  - IU^x  « a + la  (17) 


We  notice  that  this  pressure  jump  expression  has 
a sguare  root  singularity  at  the  leading  edge  of 
the  airfoil.  The  guantity,  however,  is 
Integrable  over  the  wing  chord  and  Is  used  guite 
effectively  to  yield  guantitative  estimates  of 
global  guantities  such  as  sectional  lift, 
moment,  thrust,  power  input  necessary  to  sustain 
the  motion,  energy  loss  due  to  vortex  shedding 
at  the  trailing  edge,  etc.  However,  the  result 
is  useless  as  a means  for  providing  flow  detail 
In  the  vicinity  of  the  leading  edge.  The  reason 
for  this  of  course  is  clear.  In  the 
neighborhood  of  the  leading  edge  there  occurs  a 
stagnation  point.  As  a conseguence,  the 
perturbation  does  not  remain  small  compared  with 
the  on-set  flow  as  Is  regulred  by  the  linear 
theory.  The  linear  theory  Is  therefore  not 
uniformly  valid  and  breaks  down  In  the  region 
surrounding  the  airfoil  nose.  Our  first 
objective  Is  to  correct  this  deficiency  of  the 
linear  airfoil  theory  by  determining  the 
appropriate  correction  for  the  construction  of  a 
uniformly  valid  first  order  solution  to  the 
Inviscid  unsteady  airfoil  problem. 


where  Us  and  Up  are  guantltltes  to  be  determined 
and  a,*  are  the  velocity  potential  function  and 
stream  function,  respectively. 

Substituting  Eg.  (16)  Into  Eg.  (17)  yields 


(18) 


Fig  2- Airfoil  LE  How  Fig.  3- Mapped  LE  Flow 
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From  the  above  expression  we  can  calculate  the 
nonlinear  pressure  jump  across  the  nose  of  the 
airfoil.  This  we  do  by  using  the  Bernoulli 
equation  for  unsteady  incompressible  flow 

M(q2)  + H = 9(t) 

Across  the  parabolic  airfoil  nose  the  pressure 
jump  is  Ap  = p“-  p . Therefore, 

ip/p  = - 44>t-  i(q2/2) 

The  local  time  variation  of  the  potential 
function  $ can  be  obtained  directly  from  the 
real  part  of  Eq.  (18).  This  equation  can  also 
be  used  to  obtain  q since  q = Idf/dZl. 
Consequently,  ip  across  the  nose  of  the  parabola 
V = tyJT  is 

ip/p  = 20  4UsUpVx/(l+4X)  (19) 

This  expression  is  the  leadiny  edye  counterpart 
to  Eq.  (12).  Notice  that  it  does  not  break  down 
at  the  leadiny  edye. 

MATCHING 


To  obtain  a uniformly  valid  first  order 
approximation  to  the  pressure  jump  across  the 
winy,  Eqs.  (12)  and  (19)  should  be  matched  in 
some  overlappiny  region  where  both  are  presumed 
valid.  In  Eq.  (12),  if  we  make  the  substitution 
E - 1 + x and  take  the  limit  as  E tends  to 
zero,  then 

Ap/o  ~V27e  UaQ  + 2[(c1-c2)  - UaQ/4J^2T  ( 2U ) 

In  Eq.  (19),  set  X * E and  let  E tend  to 
infinity.  This  limitiny  process  results  in 


ip/o  - - UJUp/'Vfr+  20pVT  (21) 

Compariny  coefficients  of  the  E"1^2  term  gives 
Us  » U (22) 

Ujt)  * - yJT a (t)  l 

r-  jut  / (23) 

U (t)  « V?  UU+iule3"1  ’ 

where  p 


X « oE2  + 2Fo1  - Ga2 
u * -oEj  *■  2G<Jj  + Fa2 


a2  • oEj  ♦ 2E2  ♦ oEq 

Consequently,  to  leading  order,  the 

magnitude  of  the  stagnation  point  flow  is  equal 
to  the  uniform  on-set  flow  U.  This  we  expected 
on  the  basis  of  the  steady  flow  analog  to  this 
problem.  See  Van  Dyke's  book5*,  §4.9.  An 

interesting  result  is  that  the  parallel  flow 
about  the  airfoil  leadiny  edye  is  directly 
related  to  the  strenyth  of  the  leading  edye 
suction  a0(t).  This  term  is  the  only  quantity 
that  contains  the  hist  ry  of  the  motion.  Such 
motion  history  is  due  to  vortex  sneddiny  at  the 
airfoil's  trailing  edye.  Thus  to  leadiny  order, 
the  flow  about  the  leadiny  edye  is  driven  by  the 
dynamics  and  kinematics  of  vortex  sheddiny  at 
the  trailing  edye. 


Tr  construct  a uniformly  valid  first 
approximation  for  the  pressure  jump  across  the 
airfoil  surface  we  add  the  inner  solution  (19) 
to  the  outer  solution  (12)  and  subtract  the 

common  part.  The  result  is 
£ = UaQ  [ ~=  + Vi«  S(x.t)  (25) 

where 

S(x.t)  = 2(c1«-xc. - {T  £aQ- 

Notice  that  there  is  no  singularity  in  the  above 
expression  in  the  range  |x|  < 1 where  ip  is 
evaluated. 


LAMINAR  BOUNDARY  LAYER  PROBLEM 

As  a result  of  the  asymptotic  matchiny 
technique  employed  for  this  problem,  we  have 
determined  that  the  flow  about  the  airfoil's 
leadiny  edge  oscillates  in  direct  proportion  to 
the  strenyth  of  the  leading  edye  suction.  The 
actual  flow  at  the  edye  of  the  boundary  layer 
about  the  airfoil  nose  can  be  estimated  from  Eq. 
(17)  by  taking  the  derivative  df/d<  * u.  - iu2. 
Here  uj  and  U£  are  the  velocity  components  alony 
the  (0,a)  axes  of  <t.  lie  obtain,  for  B * 0,  the 
local  flow  alony  the  parabolic  surface.  That 
is. 


Uj  * 0 

u?  * Up  + ZoUs  (26) 

Therefore,  as  one  moves  alony  the  parabola 
(either  positively  or  negatively  away  from  a 
equal  zero)  the  mean  speed  increases  in 

magnitude. 

When  the  on-set  stream  does  not  oscillate 
but  the  surface  oscillates,  the  situation 
differs  from  the  oscillatiny  dividiny  streamline 
case  only  by  the  superposition  of  a uniform, 
though  non-constant  transverse  velocity  which 
has  no  effect  on  the  relative  motion  (cf.  Ref. 
6).  Taking  advantaye  of  these  facts,  the 

relevant  boundary  layer  problem  to  consider  is 
that  of  a two-dimensional  flow  ayainst  an 
infinite  flat  plate  normal  to  the  free  stream 

where  the  plate  makes  transverse  oscillations  in 
its  own  plane.  This  is  a classical  problem  in 
boundary  layer  theory®  whose  exact  solution 
depends  on  a set  of  ordinary  differential 

equations  containing  the  reduced  frequency  o as 
a parameter.  To  conform  with  standard  notation 
for  this  problem,  we  reuse  some  notation  already 
used  for  another  purpose.  In  as  much  as  the 

principal  results  of  the  analyses  thus  far  are 

embodied  in  Eqs.  (22)  - (24)  which  are 

independent  of  coordinate  system  --  no  confusion 
will  arise. 

He  now  Introduce  a Cartesian  coordinate 
system  (x,y)  with  the  x-axis  alony  a flat  plate 
and  the  y-axis  normal  to  it  so  that  x ■ 0 is  the 
dividiny  streamline  in  the  steady  flow  outside 
the  boundary  layer  on  the  plate.  Let  (u,v)  be 
the  corresponding  velocity  components.  Outside 
the  boundary  layer  suppose 
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u ■ cx  as  y * (27) 


Let  the  plate  oscillate  along  the  x-axis  so  that 
u * ae^ut  at  y * 0 (28) 


where  c and  u are  real  constants.  The 

amplitude  of  the  plate's  speed  'a'  Is  here  a 
complex  constant.  It  is  understood  that  the 
real  part  is  to  be  taken  for  all  physical 
quantities. 


have 


Comparing  Eqs.  (22),  (23),  (26)  - ( 28)  we 


2U 


c 

a ' "R  T '“I 
The  boundary  layer  equations  are 


an  + ia.  « 


/T  U(X  + ip) 


(29) 


ut  ♦ uux  + vuy 
ux  + vy  - ° 


C X + vu 


yy 


where  u is  the  kinematic  viscosity, 
equations  are  to  be  solved  with 


aejut,  v 


cx 


0; 


These 


A similarity  solution  is  known  to  satisfy 
the  problem.  The  solution  form  is 


u • 2Uxk'  ( n)  + aejutc(n) 
v • *(2Uv)1/2k(n) 

C(n)  * cR(n)  ♦ itj(n) 
n * /Fy 


(30) 


R * 211/ v;  Reynolds  number  based  on 
Chord  length 


where  k(n),  CR(n),  c.(n)  satisf 
differential  equations.  That  is, 


ordinary 


k'"  + kk"  ♦ k'k*  ♦ 1 ■ 0 I (3i) 

k{ 0)  ♦ k'(0)  « 0,  k’(-)  - l 

kc’R  • k'cR  ♦ «,/2  • 0 \ 

Cj'*  kc'j  - k'Cj  - otR/2  • 0 > (32) 

<^(0)  » 1,  Cj(0)  ■ CR(-)  ■ Cj(-)  ■ 0 ' 


The  nonlinear  k-probl«m  is  the  classical 
Hiemenz  stagnation  point  flow.  Note  that  the 
llrear  c-problem  depends  on  the  k.solution  end 
on  the  reduced  frequency  o as  a parameter.  The 
C- solution  is  valid  for  all  values  of  o and 
'amplitude’,  a.  Another  feature  worth  noting  is 
that  the  unsteady  part  of  the  solution  is 
Independent  of  position  x.  Consequently,  the 
unsteady  part  of  the  solution  can  be  effectively 
decoupled  from  the  steady  solution. 


To  solve  the  boundary  value  problems  (31) 
and  (32)  the  differential  equations  were  written 
as  a system  of  first  order  differential 
equations  An  approximate  solution  to  the  non- 
linear k-problem  was  obtained  by  ‘Newton's 


method.  both  differential  equations  weie 
approximated  by  the  Centered-Euler  method.  The 
solutions  obtained  are  second  order  accurate. 
Figure  (4)  presents  k,  k'  and  k"  as  a function 
of  n.  Figures  (5)  - (8)  present  the  real  and 
imaginary  points  of  c and  c'  as  a function  of 
n for  select  values  of  reduced  frequency,  o. 
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LEADING  EDGE  SEPARATION 

Accordiny  to  Moore-Rott-Sears11^11  the 
point  of  separation  of  a boundary  layer  adjacent 
to  a moviny  surface  occurs  when  the  velocity  and 
the  shear  stress  simultaneously  vanish.  That 
is,  when  u * uv  * 0.  When  the  M-R-S  conditions 
are  applied  toyEq.  (30)  at  the  airfoil  nose,  x » 
0 and  the  time  dependence  is  eliminated  from  the 
resultiny  expressions,  we  obtain 

q(n;o)  ; cRCj-  cRCj  » 0 (33) 

Here,  Q is  a function  of  n that  depends 
parametrically  on  the  reduced  frequency,  o. 

For  any  specified  value  of  reduced 
frequency,  if  an  n-root  can  be  found  to  the 
equation  Q(n;o)  » 0,  then  dynamic  separation  at 
the  airfoil  nose,  x * 0,  is  believed  to  occur. 
Equation  (33)  has  been  plotted  for  a ranye  of  o. 
The  indication  is  that  dynamic  stall  does  not 
occur  at  the  airfoil  nose  for  any  value  of  o 
accordiny  o this  criterion.  See  Fiyure  (9). 
This  is  not  surprisiny  since  the  x » 0 case  is 
strictly  a shear  wave  and  symmetry  rules  out 
both  u and  uy  simultaneously  vanishing  except  at 
the  edge  of  the  boundary  layer,  n » ». 

Applyiny  the  M-R-S  conditions  when  x a 0 

gives 

2Uxf' f"  » - f"e,cos  ut  + f"a,sin  wti 

1 2 [(34) 

2Uxf'  f"  * - f*  8‘jcos  ut  ♦ f e^sin  ut  I 

where 

Bi  ’ Vr  * 4i  ci 
02  * Vi  4 Vr 

» dBj/dn,  aj  ■ dBj/dn 
Eliminatiny  x from  Eqs.  (34)  yives 
(aR<y  ajO^cos  ut  » (a^-  a^bin  ut  (3b) 
where 


“l  * V"  * CRf< 

“2  * cif'  • V" 


any  values  of  a and  n,  a time  t can  be  found 
such  that  the  M-R-S  conditions  are  met  for  the 
specified  harmonic  motion.  The  Implication  is 
that  the  M-R-S  conditions  do  not  provide  a 
useful  separation  criterion  for  this  problem. 


We  propose  the  following  methodology  for 
predicting  separation  for  this  problem: 


(i)  Set  the  steady  part  of  u equal  to  the 
amplitude  of  the  unsteady  part.  That 

is. 


2Uxk‘  (n)  ■ |a||c|  (36) 

(ii)  Set  the  steady  part  of  the  shear  stress 
(which  is  proportional  to  uv)  equal  to 
the  amplitude  of  the  unsteady  part  of 
uy.  That  is, 

2Uxk"(n)  » ja||c'  | (37) 


Eliminating  x from  Eqs.  (36),  (37) 

gives  the  leadiny  edye  dynamic  stall 
condition.  That  is, 


k‘ 1 /k‘ 


UR)2  ♦ (c1,)2 


T “7 


L K 


+ c 


I J 


1/2 


(38) 


The  right  hand  side  of  Eq.  (38)  depends 
parametrically  on  the  reduced 
frequency,  o.  For  any  specified  value 
of  o an  n-value  can  be  found  satisfying 
this  equation.  Figure  (10)  presents  a 
graph  of  the  n-root  of  Eq.  (38)  as  a 
function  of  o. 


(ill)  The  x-location  of  the  separation  point 
is  obtained  from  the  expression 

* * l4H^nrt)l/2Ull'(nrt)  (39) 

Notice  that  the  amplitude  j a | of  the 
unsteady  motion  comes  into  the  x- 
locatlr  . of  the  separation  flow. 

( i v ) Separation  occurs  for  the  motion  Oien 
the  value  of  x obtained  fro*  Eq.  (39) 
is  less  than  or  equal  to  r,  where  r is 
the  radius  of  curvature  of  the 
airfoil's  leadiny  edye.  That  is,  when 


This  equation  implies  that  for  x * 0 and  for 


|xj<  r SEPARATION 
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Fig ,iO-S«p*r*tlon  y Variation  with  { 

This  last  part  of  the  procedure  has  been 
included  since  the  local  Boundary  layer  flow  is 
not  a valid  representation  of  the  flow  adjacent 
to  the  entire  airfoil  but  only  in  a region  of 
the  order  of  the  radius  of  curvature  of  the 
airfoil  nose. 

CONCLUDING  REMARKS 

A methodology  has  been  proposed  for 
predicting  leading  edge  separation  due  to  the 
accelerating  flow  about  the  well-rounded  nose  of 
an  airfoil.  The  fluid  accelerations  are  caused 
oy  the  curvature  of  the  leading  edge  geometry 
and  the  forced  transverse  oscillations  of  the 
airfoil.  The  analyses  leading  to  an  unsteady 
separation  criterion  couples  the  gross  features 
of  the  airfoil's  transverse  motions  with  the 
details  surrounding  a laminar  boundary  layer  in 
the  vicinity  of  its  dividing  streamline. 

The  separation  criterion  has  not  been 
validated  by  comparison  with  any  experimental 
data.  This  clearly  remains  before  the  procedure 
can  be  seriously  advanced  as  a useful  tool  for 
predicting  leading  edge  dynamic  stall. 
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